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Lifshitz points are multicritical points at which a disordered phase, a homogeneous ordered phase, 
and a modulated ordered phase meet. Their bulk universality classes are described by natural gen- 
eralizations of the standard (f> 4 model. Analyzing these models systematically via modern field- 
theoretic renormalization group methods has been a long-standing challenge ever since their intro- 
duction in the middle of the 1970s. We survey the recent progress made in this direction, discussing 
results obtained via dimensionality expansions, how they compare with Monte Carlo results, and 
open problems. These advances opened the way towards systematic studies of boundary critical 
behavior at m-axial Lifshitz points. The possible boundary critical behavior depends on whether 
the surface plane is perpendicular to one of the m modulation axes or parallel to all of them. We 
show that the semi-infinite field theories representing the corresponding surface universality classes 
in these two cases of perpendicular and parallel surface orientation differ crucially in their Hamilto- 
nian's boundary terms and the implied boundary conditions, and explain recent results along with 
our current understanding of this matter. 



I. INTRODUCTION 

Lifshitz points (LP) are a particular kind of multicriti- 
cal points at which a disordered phase meets both a spa- 
tially homogeneous ordered phase as well as a modulated 
ordered one 0, 0. LI They were introduced in 1975 
by Hornreich, Luban, and Shtrikman [f|, though appar- 
ently discovered independently by two other groups 
(cf. Ref. H p. 59]). Their discovery t riggere d consider- 
able theoreticaljlM |1 E [H E El EMl EMlJll 
Il9ll20ll2ll2ll2al2^| and experimental interest l25ll2rf. 
which has continued over the years, and after a phase of 
somewhat reduced intensity, has regained a lot of momen- 
tum recently [H l28l [29L |30l |3lL 132. 133 — in particular, 
on the theory sidejMM]^!!! THU El EJ El El 

El El El EE El IMlanH^r 

The physics of LP is interesting for a variety of reasons. 
Let me mention a few. 

(i) A wealth of physically distinct systems exist that 
are either known to have LP or for which LP have 
been discussed; this includes such diverse systems 
as magne ts |25l. ferroelectrics [31], polymer mix- 
tures |28l l38l l52| , liquid crystals [53j , systems un- 
dergoing structural phase transitions or domain 
wall instabilities [54| , organic crystals d US > an d 
even superconductors |56j . 

(ii) The physics of LP embodies many of the crucial 
concepts of the modern theory of phase transitions 
and critical phenomena, yet has been explored to a 
much lesser degree than critical behavior at con- 
ventional critical points. The best studied uni- 
versality classes of bulk critical behavior are the 
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ones for d-dimensional systems with short-range 
interactions and an rt-component order parame- 
ter field cf), represented by the 0(n) 4> 4 models. 
For them, very detailed — and in part impressively 
accurate — results have been worked out by means 
of sophisticated renormalization group (RG) ap- 
proaches [57l IHfij , series expansions 1 and com- 
puter simulations |60| ; and many of these theoreti- 
cal predictions have been checked by careful exper- 
iments. 

By contrast, the application of modern field- 
theoretic RG approaches to the study of criti- 
cal behavior at LP is a fairly recent development 
IH El El EH The two-loop RG analysis of criti- 
cal behavior at m-axial LP in d = — e dimen- 
sional systems Shpot and myself El 03 managed 
to perform for general values < m < d has finally 
yielded the e expansions of all critical exponents to 
second order. The estimates obtained by means of 
these series expansions for the values of the critical 
exponents for the scalar uniaxial case n = m = 1 in 
d = 3 dimensions agr ee q uite well with up-to-date 
Monte Carlo results [61|. Unfortunately, we are 
aware of only a few high-temperature series esti- 
mates [ll |2l |24j , none of which is very recent. On 
the experimental side, renewed activity is notice- 
able. Aside from the recent work on polymer mix- 
tures jH new experiments on magnetic sys- 
tems have been reported [2l|3(|. However, so far 
the latter have not produced results for the critical 
and crossover exponents of the m = n = 1 LP point 
of significantly greater accuracy than achieved in 
previous studies |HE3- 

(iii) Compared with critical points, LP provide ad- 
ditional challenges. Since they are multicritical 
points, a further thermodynamic variable besides 
temperature T must be fine-tuned to reach them. 
Furthermore, precise experimental investigations of 
their critical behavior should include verification of 



2 



the expected crossover scaling forms and expected 
to involve the choice of proper nonlinear scaling 
fields 

On the theoretical side, progress in analytical RG 
analyses has been hampered by the substantial 
technical difficulties one encounters in computa- 
tions of Feynman diagrams beyond one-loop or- 
der. The origins of these problems are twofold: 
the anisotropic nature of scale invariance that holds 
at the LP, which implies that the free propagator 
does not reduce to a simple power at the LP but 
involves a scaling function; and the fact that this 
scaling function in position space turns out to have 
a rather complicated form in g eneral |43j . The 
progress made recently 0, IHF m handling such 
field theories could pave the way for systematic in- 
vestigations of general aspects of anisotropic scale 
invariance (ASI) in systems with short-range in- 
teractions. One imp ortant question that has been 
raised long ago jHS 153 but not yet answered in 
a truly convincing fashion is the following. Scale 
invariance, in conjunction with translation and ro- 
tation invariance, and short range of interactions, is 
known to normally imply invariance under a larger 
symmetry group, n amely under conformal trans- 
formations |65t l66t Wfy . Does ASI likewise entail 
the invariance under additional nontrivial continu- 
ous transformations? Henkel has played with this 
idea for years making concrete proposi- 

tions for transformations under which two-point 
correlation functions should be invariant, he has 
come up with definite predictions for the form of 
the associated scaling functions, which appear to 
be consistent with Monte Carlo results [61j for 
the three-dimensional ANNNI model, yet remain 
to be carefully checked by analytical calculations 
[6^. FFof . The field theories representing the univer- 
sality classes of critical behavior at m-axial LP are 
particularly well suited for such scrutiny, not least 
because the parameter m can be varied. 

(iv) Since LP involve both modulated ordered phases as 



well as ASI, rich and interesting boundary critical 
phenomena 0, Iz3| mav be expected to occur 
near them. The systematic investigation of such 
phenomena, in particular, via field-theoretic RG 
tools, is still in its infancy [63. t74l r75l 17a . 1771 r7^| . 

In this contribution, I will briefly survey the progress 
made recently in the application of field-theoretic RG 
methods to bulk and boundary critical phenomena at LP, 
compare its results with those from other sources such as 
Monte Carlo simulations, highlight some of the central 
issues and difficulties, and indicate directions for further 
research. We begin in the next section by specifying the 
models, then deal with their bulk critical behavior, be- 
fore we turn in Sec. lllll to the issue of boundary critical 
behavior. 

II. CONTINUUM MODELS AND BULK 
CRITICAL BEHAVIOR 

A. Continuum models 

Having in mind systems whose microscopic interac- 
tions are either short ranged or of a long-range kind 
that is irrelevant in the RG sense, we consider contin- 
uum models with a Hamiltonian of the form 

U = [ C h (x)dV+ [ d(x)dA, (I) 

where Ch{x) and C\{x) are functions of the n-component 
order parameter 4>(x) — {(f> a (x)) and its derivatives with 
respect to the coordinates (x a ,xp) = x. We index the 
first m Cartesian coordinates by a; they refer to the Tri- 
dimensional subspace to which the modulation of order is 
confined. The remaining m = d~m ones are labeled by (3. 
When we deal with boundary critical behavior, the vol- 
ume and surface integrals dV and dA extend over 
the half-space = {x = {r,z)\r G K d_1 , < z < oo} 
and the z — hyperplane 23, respectively. To investi- 
gate bulk critical behavior, we may as well take 03 = M. d 
and forget about the boundary piece in Eq. QJ, choosing 
appropriate (periodic) boundary conditions. 



Unless stated otherwise, the bulk density is 

4,(*) = f(£W +\ E (« 2 + f £(0«*) 2 + 5* 2 + £l*l 4 (2) 

^q=1 ' p=m + l a=l 

in the sequel. Here d a = d/dx a and dp = d/dxp, and a > as well as u > is assumed. For the time being we focus 
on bulk critical behavior. Let us therefore postpone the choice of the boundary density L\ to Sec. Mil Our selection 
P)l of Ch reflects two tacitly assumed properties: 0(n) invariance, and isotropy in the m-dimensional a-subspace of 
coordinates. An investigation of the effects of spin anisotropies breaking the O(n) invariance of Ch may be found in 
Ref. |l8j |: they will not be considered here. However, the role of "space anisotropies" reducing the Euclidean invariance 
in the a-subspace [HlJ will be briefly discussed at the end of this section. 



From Eq. it is easy to understand how a LP can on T and a second thermodynamic variable, a non- 
occur. The interactions constants <x, . . . , u all depend 
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ordering field g Q such as pressure (charge-transfer salts 
|sll55|L a ratio of next-nearest neighbor (nnn) antiferro- 
magnetic and nearest-neighbor (nn) ferromagnetic inter- 
actions along an axis (ANNNI model 0), or a magnetic 
field component in the subspace orthogonal to the or- 
der parameter (the orthorhombic magnetic crystal MnP 
US lH IM E2). Assuming that the coefficient of the 
{dp4>) 2 term does not change sign, we have absorbed 
it in the amplitude of <fi so that it becomes 1/2. Lan- 
dau theory gives a disordered phase for f > provided 
p > 0, separated from a homogeneous ordered one by 
the critical line r c (p > 0) = 0. For negative p, a con- 
tinuous transition from the disordered to a modulated 
ordered phase occurs across the so-called "helicoidal sec- 
tion" t = f c (p < 0) of the critical line, which joins the 
"ferromagnetic section" at the LP f — p — (see, e.g. 
Fig. 1 of Ref. Q). The other phase boundary emerging 
from the LP separates the homogeneous ordered from the 
modulated ordered phase. The transitions across it can 
be of first or second order; for cases with a scalar order 
parameter they are generically discontinuous, whereas for 
specific models with a vector order parameter they turn 
out to be continuous 1791. 



B. Critical exponents, anisotropic scale invariance 

In Landau theory, the helicoidal section fh c = f c (p < 
0) varies as fhc ~ P 2 near the LP. Beyond Landau theory, 
the LP and the phase boundaries — supposing they still 
exist — get shifted as a result of fluctuations, and the he- 
licoidal section of the critical line is expected to behave 
near the LP as 

fhc (Sp) - f L p = Sf hc ~ l^l 1 ^ ~ \8g\ x ^ . (3) 

Here Sp and Sg denote deviations of p and g from their 
values plp and g^p at the LP. We have introduced the 
crossover exponent ip, whose mean- field value is <pmf = 
1/2, and utilized the fact that 5p ~ 5g near the LP. 

In the modulated ordered phase, the order is modu- 
lated with a wave vector q mo d (T, g) depending on T and 
g. Since homogeneous order corresponds to q mo d = 0, 
Qmod must also vanish at the LP. Its limiting behavior as 
the LP is approached along the critical line's helicoidal 
section T = Th c (<?) is governed by the wave- vector expo- 
nent P q , defined via 



analogs of the usual correlation exponent rj by 



9mod(T hC)ff )~ \6pf«~ \Sgf" 



(4) 



Other important critical exponents characterize the 
scale invariance at the LP. Let us set f = p = u = 
in Eq. @ and transform to momentum (q) space to 

( 2) 

obtain the two-point bulk vertex function (q) in 
the Ornstein-Zernicke approximation. Its q-dependence 
reads a (q a q a ) 2 + qpqp, where repeated indices a and (3 
are to be summed over 1, . . . , m and m + 1, . . . , d, respec- 
tively. Beyond this classical approximation one antici- 
pates again nontrivial power laws. Hence one introduces 



2 for q a 



0. 



4 for qp = 0. 



(5) 



These relations mean that q a scales as q a ~ [qp) B , with 
the "anisotropy exponent" 



(2-77 L2 )/(4-77 L 4) ; 



(6) 

1/2. Likewise 



in Landau theory it takes the value #mf 
x a ~ (xp) 6 . 

To formulate ASI in position space, let us consider a 
perturbation go J <a O(x) dV of the fixed-point Hamilto- 
nian associated with the LP, where 0{x) is a scaling 
operator with scaling dimension A[C]. Let yo be the 
RG eigenexponent of the associated scaling field go, so 
that go — » go{ty = £~ vo go under scale transformations 
x p xpl. Since the scaling dimension A[C] and the 
eigenexponent yo must add up to minus the scaling di- 
mension of the volume V — dV, which is m + m 0, we 
have 

yo = m + m6- A[0] . (7) 
The operators 0(x) satisfy 

O(£ e x a ,tx ) =r Alo] O(x a ,x0) (8) 
(ASI) in the long-scale limit I — > 0. 

C. Field theory and e expansion 

For a conventional critical point it is known that below 
the upper critical dimension d* — 4, where hyperscal- 
ing is valid, two independent critical exponents exist in 
terms of which all critical indices characterizing the lead- 
ing infrared singularities can be expressed. They derive 
from the scaling dimensions of 4> and the energy density 
4> 2 , or equivalently, the RG eigenexponents yh and y T . 
Furthermore, there are just two metric factors, one as- 
sociated with each of the corresponding scaling fields h 
and t ("two-scale factor universality" |S3)- -"- n the case 
of m-axial LP, the upper critical dimension is |j| 



d*{ m ) = A+- . 



(9) 



The easiest way to see this is to determine the dimension 
d = d*(m) below which the Gaussian scaling dimension 
of u becomes positive; the Ginzburg criterion yields the 
same result. 

In view of the different scaling of x a and xp, and the 
need to fine-tune an additional variable — p or g — , it is 
natural to expect that four critical exponents will be re- 
quired to express the bulk critical exponents of the LP for 
d < d*{m). Of course, some of these might turn out to 
be trivial, taking on values independent of d and m. For 
example, one might anticipate the anisotropy exponent 
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9 to retain its mean- field value for d < d*(m). However, 
the e-expansion results of Shpot and myself 0, EH have 
revealed that nontrivial m-dependent contributions to 9 
appear at order e 2 . The bulk operators 0(x) from whose 
scaling dimensions the four independent bulk exponents 
derive are given in Tabled along with the associated scal- 
ing fields and their RG eigenexponents. Each of these 
four scaling fields involves a nonuniversal metric factor. 
Hence a four-scale- factor universality applies. 



Given a line of upper critical dimensions d*(m), one 
should be able to expand about any point on it. Al- 
though this goal was identified at a very early stage p|, 
its implementation turned out to be very demanding and 
took a long time. In Refs. |43|,|46j a two-loop RG analysis 
was performed in d* (to) — e dimensions for general values 
of to. This gave the e expansions of the four independent 
bulk critical exponents t]l2, 0, vl2, and tp, as well as the 
correction-to-scaling exponent oj u , to 0(e 2 ). 



Technically, a massless minimal-subtraction renormalization scheme was employed. To define the ultraviolet (uv) 
finite renormalized theory, the reparametrizations 

= z l /2 ^rcn , cr = Z a a , u a~ m/A F m<e = p e Z u u , 
t — tlp = p 2 Z T [t + A T p 2 ~\ , (p-Plp) ct~ 1/2 = pZ p p . (10) 

were made, where p is a momentum scale, while F m ^ is a convenient (uv finite) normalization factor whose precise 
choice need not worry us here. All renormalization factors Z^, Z G , Z T , Z p , and Z u were computed to 0{u 2 ). From the 
result for Z u the RG beta function (3 u (u, e) follows to order it 3 ; the other Z-factors yield RG functions whose values at 
the nontrivial root u*(m, e) of (3 U give the critical exponents. The main consequence of the renormalization function 
A T is that the scaling field with the RG eigenexponent 1 / vl2 becomes a linear combination of r and p 2 0, |6^| • 

What makes calculations beyond one- loop order complicated is that the scaling function $ m ,d{v) of the free bulk 
propagator at the LP, 



G b (\x\ 



d d q exp(iq • x) 



(2ir) d / 2 qpqp + a{q a q a ) 2 

m/4 (xfixfi)- 1 ^' 2 $ m>d (v) , v = {bxpXfiY 11 ^^-^ , (11) 



is a difference of generalized hypergeometric functions. While these increase in general exponentially as v — ► oo, their 
difference has an asymptotic expansion in inverse powers of v that does not terminate except for special choices of 
(to, d), such as (2,5) and (6,7), where it reduces to elementary functions 0,0- Therefore, the two-loop series 
coefficients of the renormalization factors could not be computed analytically for general to. However, they — as well 
as the implied e-expansion coefficients of the critical exponents — could be written in terms of four single integrals 
j<j,{m), Jo-(to), j p (m), and J u (m) of the form J Q dv f(v;m), where f(v;m) involves ^ m ,d'{m){ v )i analogous (related) 
scaling functions, and powers of v ^(J- For to = 0, 2, 6, 8, these integrals could be computed analytically; for other 
values of to they had to be determined by numerical means. 



The resulting e expansions of the critical exponents 
A = vl2, ■ ■ ■ ! P and the correction-to-scaling exponent 
lu u take the form 

A(n, to, d) = A MF + Ai(»)e+ \ 2 {n,m)e 2 + 0(e 3 ) . (12) 

Note that Ai(n) is independent of to, so that the in- 
dependence starts at order e 2 . This means that the co- 
efficients Ai(n) coincide with their to = counterparts 
for the standard </> 4 model for all exponents that remain 
meaningful when m — 0. (Recall that exponents such as 
?7L4, tp, and 9 are not needed in the isotropic case to = 0.) 

The result (I12|) allows several interesting checks. First, 
if we substitute the analytically known m — values of 
the integrals j L (m) into it, choosing A = «L2, ^L2, and 
oj u , then the familiar expansions to 0(e 2 ) of the standard 
exponents r], v, and uj u of the ip 4 model are recovered. 
A second check concerns the special cases to = 2 and 
to = 6. Owing to enormous simplifications, the two-loop 



RG analysis can be performed fully analytically. The re- 
sults one obtains in this fashion are fully consistent with 
what one gets upon insertion of the analytically known 
values of j L (2) and j L (6) into the two- loop expressions for 
general to. Third, considering the case of the isotropic 

lp nm, 

one can set d = to = 8 — eg in Eq. i|12[l and 
expand to second order in eg = 2e. The limiting values 
it (8—) are again known analytically 0,^^]. Considering 
exponents that remain meaningful in the isotropic case 
to = d, such as 77^4 or 1/^4 = 9 vl2, we can derive their ex- 
pansions in eg to 0(e% 2 ) from Eq. (|12fl . The results agree 
with those obtained via a direct analysis of the isotropic 
model with to = d in 8 — eg dimensions 0, |48|| . 

A cautionary remark is appropriate here. As a can- 
didate for an experimental system with an isotropic Lif- 
shitz point, ternary mixtures of A and B homopolymers 
and AB diblock copolymers have been studied both ex- 
perimentally [23] and theoretically jE^]. In their case, 



5 



TABLE I: Bulk scaling operators O(x), associated scaling dimensions A[0], bulk scaling fields go, and their RG eigenexponents 
yo, giving the four independent bulk critical exponents of the LP. 

Q(x) AJO] go yo 

4> {m + me -2 + J? i2 )/2 h (m + md + 2 - T] L2 ) /2 

{d a d a (f>) 2 fh + md- 46 + 2 a 46-2 

4> 2 m + m6 — 1/ul2 = 0- — <*l)/vl2 T 1/VL2 

(d a (j))d a 4> fh + md — ip / i/L2 p <f/vL2 



modulated order occurs in the lamellar phase. While 
self-consistent field theory predicts the transition from 
the disordered to the lamellar phase to be continuous 
|52j . theoretical arguments in favor of a first-order tran- 
sition have been presented |8l|. This would mean that 
there is actually no isotropic LP. According to some ex- 
periments (see the discussion in Sec. 7 of Ref. |5^| ). the 
Lifshitz point found in mean- field (MF) theories gets ap- 
parently destroyed. Unfortunately, recent Monte Carlo 
simulations |52| were not able to decide whether the tran- 
sition between the disordered and lamellar phases is of 
first order or continuous. However, they yielded modi- 
fications of the MF phase diagram similar to those seen 
in experiments — in particular, no LP. If fluctuations in- 
deed preclude the appearance of an isotropic LP, then 
the e expansions of the critical exponents for that case 
are mainly of academic interest. Nevertheless, their con- 



sistency with the results for general m is very gratifying 
from a mathematical point of view. 

The series-expansion results for general m can be, and 
were, used in particular to obtain approximate values for 
the critical exponents of the uniaxial LP with n = 1 at 
d = 3 S3| . Both experiments on MnP [25| well as Monte 
Carlo calculations for the ANNNI |20L l61| model provide 
clear evidence for the existence of such a LP. Recent field- 
theoretic estimates are v^i — 0.75, (3l = VL2 A[</>] ~ 0.22, 
9 = v Li /v L2 ~ 0.47, ip ~ 0.68, a L ~ 0.16, and 7L ~ 1.4 
|46j . The agreement with current Monte Carlo results, 
which gave a L = 0.18 ± 0.03, j3 L = 0.235 ± 0.005, and 
7l = 1.36 ± 0.03, is fairly good. For more detailed com- 
parisons covering also other cases, experimental work, 
and further theoretical estimates the reader is referred to 
Refs. HEllll. 
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D. Space anisotropies 

A natural generalization of the ANNNI model is the biaxial nnn Ising (BNNNI) model, which has competing nn 
and nnn interactions along two cubic axes rather than along a single one. In d dimensions, even m-axial variants 
of the latter, "mNNNI models" with m < d, can be considered. The continuum models onto which they map upon 
coarse graining generically have fourth-order derivative terms breaking isotropy in the a-subspace. Their symmetry 
may be cubic or — if we consider similar microscopic systems involving other crystal lattices — even weaker. Hence, 
whenever m > 1, the bulk density J2J) should be supplemented by anisotropic contributions of the form 



a o IIL 

4") = ^mTi% 2a3ai (d ai d a2 <p)d a3 d ai <f> = |u>$>^) 2 



(13) 



where all tensors 7 , a 1 Q 2Q3a4 permitted by symmetry must be included. The w% are dimensionless interaction constants. 
For cubic symmetry, only the first term on the far right remains. 



The effects of such space anisotropies were investigated 
in Ref. [HI]]. A new renormalization factor Z Wi is re- 
quired for each independent anisotropy, and both these 
as well as the previously introduced renormalization func- 
tions [Eq. nnj] now depend on u and the renormalized 
anisotropies mj. Specifically, the crossover exponent if2 
associated with the cubic anisotropy w was computed 
to 0(e 2 )- For m = 2 and m = 6, the 0(e 2 ) coefficient 
could be determined analytically, for other values of m 
expressed in terms of another (numerically computable) 
single integral. It turned out to be small, but positive. 



For example, for m — 2 and n = 1 its evaluation at 
e = 2 yielded the d = 3 estimate if 2 = 1/81 ~ 0.012. 
Thus the isotropic fixed point is u>i — unstable, at least 
for small e. Whenever such anisotropy is present, the 
previously found universality classes should not apply. 
Unfortunately, no new stable fixed point could be found. 
A detailed clarification of the behavior for Wi 7^ re- 
mains a challenge. It would be interesting to investigate 
the role of such anisotropies in Monte Carlo simulations 
of suitably designed three-dimensional models (e.g., the 
BNNNI model), albeit deviations from the Wi = uni- 
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versality classes may be difficult to measure because of 
the smallness of (fi2 ■ 



III. BOUNDARY CRITICAL BEHAVIOR AT LP 

The study of boundary critical behavior at LP started 
with Gumbs' work based on Landau theory |74| . in which 
z was taken to be an a-direction. Later considerably 
more detailed MF analyses [TBI and Monte Carlo 
calculations [77| of semi-infinite ANNNI models with 
perpendicular {z = a-direction) and parallel (z = [3- 
direction) surface orientations were performed. So far, 
detailed field-theoretic RG studies were made only for 
the case of parallel surface orientation [(J^, ■ 

Let me emphasize that the two primary types of surface 
orientations (|| or _L) correspond to substantially distinct 
cases. This can be seen from the following observations: 
First, z scales differently, namely, as £~ x and £~ e , respec- 
tively. This has an immediate consequence. Consider a 
perturbation g^m J^G^r) dA, where s (T*) is a bound- 
ary operator with scaling dimension A^* 8 ] and hence 
has the ASI property JSJl . The analogs of Eq. Q for the 
RG eigenexponent vqh of gQ<s differ depending on the 
surface orientation: 



m + m 9 — A 



II.-U 



A[0 



<8i 



A i 



e. 



(14) 



Second, owing to the different engineering dimensions 
[z] = and [z] — ct 1 / 4 /! -1 / 2 , power counting consid- 
erations to estimate the relevance or irrelevance of con- 
tributions to the surface density L\ differ. Third, since 
Ginzburg-Landau theory yields differential equations for 
the order parameter of second (||) or fourth (_L) order in 
d z , either a single or else two boundary conditions are 
needed at z = and z — oo. 

To bring the problem into focus, let me recall that in 
the m = case of the standard semi-infinite <fi model it is 
sufficient to choose L\ = \c(j) 2 , unless terms breaking the 
0(n) symmetry are permitted [Z2 (which will b e avoided 
here) . On the basis of power counting alone, one might 
think that the symmetry-allowed monomial (pd n <p (where 
d n means derivative along the inner normal), should be 
included as well. But this is redundant because of the 
boundary condition d n <p = ccp, which as usual follows 
from the boundary part of the classical equation SH = 
and holds beyond Landau theory inside of averages. 

The surface enhancement variable c determines the 
type of surface transition that occurs at bulk criticality: 
Depending on whether its deviation Sc = c — c sp from a 
special value c sp satisfies Sc > 0, Sc = or Sc < Oan ordi- 
nary, special or extraordinary transition occurs [7H l72j | , 
provided the dimension of the surface, d — 1, exceeds the 
value below which a long-range ordered surface phase in 
the presence of a disordered bulk is not possible (i.e., if 
d > 2 and d > 3 in the Ising and n > 1 cases, respec- 
tively). 



What modification occur in the m > LP case? They 
are easy to understand if the surface orientation is par- 
allel: An additional derivative term must be included in 
£%, which thus becomes |62| 



o i III 

a=l 



(15) 



Since [c] = [ct 1 / 2 ^] = /i, the variable Act 1 I 2 is di- 
mensionless. The implied boundary condition reads 
(d n — Xd a d a )4> = c(j>; it can be employed to con- 
clude that contributions to C\ of the form drd n <f> and 
(d a <fi)d n d a <fi are redundant. By contrast, the inclusion 
of the term oc A is necessary: Not only is it required 
to absorb uv singularities of the theory, but it would be 
generated under the RG if originally absent. This can 
be seen as follows: In order to renormalize the model de- 
fined by Eqs. QJ, and l|15l) . we must complement the 
reparametrizations (|10|) by 



Act- 1 / 2 

C Car, 



X + P x {u,X 7 e) , 

H Z c [c + A c (u, X, e) p] 



(16) 



Here the surface renormalization factors Z\ and Z c de- 
pend on u and A, just as P\ and A c . At 0(u 2 ), P\ does 
not vanish for A = 0, so a nonzero A gets indeed gen- 
erated. Furthermore, there are no RG fixed points at 
A = on the hyperplane u — u* (see Fig. 2 of Ref. [62jp. 
The fixed points associated with the ordinary, special, 
and extraordinary transitions turn out to be located at a 



nontrivial A- value A3 



Ao(m) + 0(e) and c = c* rd = oo, 



c* p = 0, and c* x = — oo, respectively. 

Before continuing our account of the available re- 
sults for this parallel case, let us briefly discuss how to 
choose C\ when the surface orientation is perpendicu- 
lar. Clearly, the two monomials included in Eq. I|15|) 
should be expected here as well, although different cou- 
plings ought to be associated with (d a <p) 2 for a = 1 (z- 
direction) and a > 2. As long as terms breaking the 0(n) 
symmetry can be ruled out, the choice 



L 1 - 



x m \ 
i 2 + ^(a^) 2 + &H0+y(^^ (17) 

a=2 



should be sufficient. From the vanishing of the contri- 
butions J(g • • • Sd n <p and . . . S<p to SH. two boundary 
conditions on 03 are found, namely 



adl + (b 



p)d n +c^-X^dl 



= 0, 



-<t^ + A_l+6 = 0. (18) 



They tell us that the monomials cf)d 2 <p, (d n <fi)d 2 4>, and 
0(9^ </> (which are potentially dangerous for e > accord- 
ing to power counting) are redundant. A detailed RG 
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analysis of the model with the bulk and surface densities 
© and l(T7|) remains to be done. 

In the case of parallel surface orientation, it is pos- 
sible to investigate the ordinary transition without re- 
taining the dependence on A and c 62]: In the limit 
c — > c* rd = oo a Dirichlet boundary condition applies 
and the dependence on A drops out (resides only in met- 
ric factors). Hence one can set c = oo and A = 0, choos- 
ing from the outset Dirichlet boundary conditions for the 
bare theory. The critical exponent [3\ of the surface or- 
der parameter s (r) = 0(r,O) follows via the boundary 
operator expansion 

<Kr, z) « C(z) d n cj> , C{z) ~ z ^]-^] 5 (19) 

z— >u 

giving (3° /vl2 = A[d n 4>]. Hence one must study multi- 
point cumulants involving an arbitrary numbers of fields 
4> and boundary operators d n cf>. This strategy was fol- 
lowed in Ref. |62j and utilized to determine the critical 
index /3° rd to 0(e 2 ) for general < m < 6. The e 2 term 
involves a further single integral j\ (m) , which again could 
be computed analytically for m = 0, 2, 6, though only nu- 
merically for other values. All other surface exponents of 
the ordinary transition can be expressed in terms of a 
single one, e.g. /3° rd and four independent bulk indices. 
The form l|12fl of the e expansion, with m-independent 
0(e) terms, also applies to these surface exponents. Fur- 
thermore, for m — > their expansions to 0(e 2 ) turn into 
the known ones [T^, of the standard semi-infinite cj> 4 
model. The d = 3 estimates one obtains from these e 
expansions in the uniaxial scalar case m = n = 1 (e.g., 
/3° rd ~ 0.68... 0.7) agree reasonably well with recent 



Monte Carlo results for the ANNNI model [73, which 
gave /3f d = 0.687(5). 

The special transition is harder to analyze because the 
A-dependence must be retained, though c can be set to 
its fixed-point value c* p = 0. A recent one-loop analysis 
|78j showed that 0^ agrees with the bulk exponent (3l to 
O(e) and that the crossover exponent $ associated with 
c becomes m-dependent alr eady at O(e). According to 
recent Monte Carlo results Ullnl Pl P = 0.23(1) and 
(3 L = 0.238 ± 0.005. Thus the difference 0i - (3 L seems 
to be small indeed. 

Returning briefly to the case of perpendicular sur- 
face orientation, let me conclude with a — hopefully 
educated — guess concerning the ordinary transition. I 
expect that the asymptotic behavior at this transition 
is described by a theory that obeys the boundary con- 
ditions cfy 8 = d n cj) = 0. The critical exponent /3\ in this 
case should follow from the boundary operator expansion 

z— >0 

(20) 

and be given by (3° vd /vl2 = A[<9 2 (/(>]. 
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